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O ', Abstract 
^ ' 

, We present a comprehensive overview of theory-data comparisons for inclusive jet pro- 

duction. Theory predictions are derived for recent parton distribution functions and 
Q>^ , compared with jet data from different hadron- induced processes at various center-of- 

' mass energies ^/s- The comparisons are presented as a function of jet transverse mo- 

mentum pt or, ahernatively, of the scahng variable xt = IpT I \fs- 

:t . 1 Introduction 

X 

' In this article we provide a comprehensive overview of theory-data comparisons for inclusive 

jet cross section measurements made in hadron-induced processes between 1993 and today. 
The comparisons are made for the central region in (pseudo-) rapidity. Data sets include 
inclusive jet cross sections measured in hadron-hadron collisions at center-of-mass energies 
\/s between 200 GeV and 7TeV and in deeply inelastic scattering (DIS) at -^/i = 300 GeV 
and 318 GeV. The theory calculations are computed consistently for all data sets using recent 
parton distribution functions (PDFs) and common choices for the strong coupling constant 
as{Mz) and the renormalization and factorization scales. Ratios of data and theory are 
presented as a function of jet transverse momenturr0 pT or, alternatively, as a function of 
the scaling variable xt = 2pT / \/s. 



^Some of the measurements discussed in this article have been made as a function of jet transverse 
energy Et- Throughout this article we use pr to refer to both, either the jet transverse momentum or the 
jet transverse energy, as appropriate. In hadron-hadron collisions, "transverse" refers to the momentum 
or energy component perpendicular to the hadron-hadron beam axis. In DIS, "transverse" refers to the 
momentum component perpendicular to the photon-proton axis (which is the z-axis in the Breit frame). 
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This article is intended as a persistent repository for the newest theory-data comparisons 
possible in jet production in hadron-induced processes and will be updated whenever new 
relevant data and/or theory results are published. 

In section 2 we discuss how the theory results are obtained while details of the data sets 
are described in section 3. Comparisons of theory and data are presented in section 4. 

2 Theory Predictions 

Conventionally, experimental measurements are corrected for instrumental effects and the 
results are provided at the "particle level" [I]. Corresponding theory predictions for jet 
observables Otheory are therefore obtained as the product of a perturbative QCD (pQCD) 
result, OpQCD, multiplied by a correction factor for non-perturbative effects, Cnon-pert 

Otheory = OpQCD ' t-non— pert ■ (1) 

For inclusive jet cross sections in hadron-hadron collision and in DIS, the fixed-order pertur- 
bative expansion is known to next-to-leading order (NLO) in the strong coupling constant 
as, which is of order 0{a1) for DIS and 0{af) for hadron-hadron collisions. For inclusive 
jet cross sections in hadron-hadron collisions, additional 0{a^) corrections from threshold 
corrections have been computed [2]. These results are used for the calculations in this arti- 
cle. Renormalization and factorization scales, fir and fj,f, are set to fir = fJ-f = Pt, where px 
is the transverse momentum (or, where appropriate, transverse energy) of the individual jet. 
In all cases fir and fif are larger than the b-quark mass, and thus all calculations are made 
for five active flavors {uf = 5). The proton PDFs are taken from the PDF parametrizations 
MSTW2008 [31, CTIO [1], NNPDF2.1 [5], or HERAPDF1.5 [HE]. Ah of these parametriza- 
tions are available in NLO and NNLO approximations (except for CTIO which are only 
available in NLO) as well as for a series of as{Mz) values. In each case, we use those PDF 
sets which were obtained for as{Mz) = 0.118 which is closest to the current world average 
value [8]. The same value of as{Mz) is also employed in the matrix elements. The NLO 
PDF parametrizations are used for the comparisons with DIS data, together with the NLO 
matrix elements. The NNLO versions (where available) are used for the comparisons with 
data from hadron-hadron collisions, together with NLO matrix elements plus ©(af) con- 
tributions from threshold corrections. Correspondingly, in these comparisons Og is evolved 
using the numerical solutions of the 2-loop and 3-loop approximations of the renormalization 
group equation for DIS and hadron-hadron collisions, respectively. All pQCD contributions 
are computed within the framework of fastNLO [9] using NLOJET++|101 111] for the LO 
and NLO contributions and code from the authors of Ref. [2] for the threshold corrections. 
The anti-ZcT jet algorithm is taken from the implementation in Fast Jet [12\. 

The factor Cnon-pert includes non-perturbative corrections due to hadronization and the 
underlying event. These corrections are usually obtained using Monte Carlo event generators 
and the values are published together with the measurement results. Whenever this was 
not the case, we have computed the correction factors ourselves using pythia 6.4 |13j with 
tune A [n]. For DIS data at high Q^, underlying event corrections are assumed to be 
negligible. In those cases, the pQCD results are only corrected for hadronization effects. 

For some of the data sets, additional theoretical effects have already been folded into 
the published results. This is taken into account as follows: 
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• For some pp data sets (as mentioned in section 3), corrections for the underlying event 
have already been applied to the data. In the corresponding calculations, pQCD 
theory is only corrected for hadronization effects. 

• Observables in DIS are usually corrected for higher order QED effects which may 
either include corrections for the running of the QED coupling a, or not. The former 
is the standard choice by the ZEUS Collaboration, and therefore the corresponding 
pQCD calculations are made for a fixed value of a = 1/137. Jet measurement results 
from the HI Collaboration are not corrected for the running of a, and therefore the 
corresponding pQCD calculations include a running value of a{Q). 

3 Data Sets 

In this section, we give a brief overview on the data sets which are used in the theory-data 
comparisons. The phase space requirements are given as stated in the cited publications 
which include further details regarding the exact definition of the variables. The recombi- 
nation schemes used in the jet algorithms are usually the E'-scheme [15] . or the E'j'-scheme 
(frequently referred to as "Snowmass convention") |16j. In some cases other recombination 
schemes have been used. The following measurements of the inclusive jet cross section as a 
function of pT are used in the comparisons. All data are corrected to the particle level as 
defined in Ref. pQ, except where mentioned otherwise. 

1. Proton-proton scattering at ^/s = 7TeV at the LHC 

(a) CMS Collaboration (2011) d^t = 34 pb^^ [17J 

Measurement of d?a/dpTdy for 18 < pt < 1100 GeV in six rapidity regions 
\y\ < 3.0 using the anti-fcy jet algorithm with R = 0.5 in the i?-scheme. In our 
comparisons we display the results for |y| < 0.5. Non-perturbative corrections are 
determined as the average of pythia 6.422 with tune D6T [18j and herwig++ 
2.4.2 [19] with the default tune of herwig++ 2.3 (as detailed in Ref. [20]). 

(b) ATLAS Collaboration (2010) d^t = 17 nb^^ [2T] 

Measurement of cPa/dpTdy for 60 < pr < 500 GeV in five rapidity regions 
\y\ < 2.8 using the anti-/cT jet algorithm with R = 0.4 in the i?-scheme. In 
our comparisons we display the results for < \y\ < 0.3. Non-perturbative 
corrections are determined using pythia 6 with tune MC09. 

(c) ATLAS Collaboration (2010) dnt = 17 nb^^ [H] 
The same as in 1. (b), but for R = 0.6. 

2. Proton-anti-proton scattering at ^/s = 1.96 TeV at the Fermilab Tevatron collider 

(a) CDF Collaboration (2008) dnt = 1.13 fb"! [22] 

Measurement of d'^a/dpxdy for 62 < < 700 GeV in five rapidity regions \y\ < 
2.1 using the Run II iterative midpoint cone jet algorithm (CDF version) with 
R = 0.7 in the E-scheme. In our comparisons we display the results for 0.1 < 
\y\ < 0.7. Non-perturbative corrections are determined using pythia with tune 
A. 

(b) CDF Collaboration (2007) dnt = LOfb^i [23] 

Measurement of d^a/dpxdy for 52 < px < 700 GeV in five rapidity regions \y\ < 



3 



2.1 using the inclusive jet algorithm with R = 0.7 in the E'-scheme. In 
our comparisons we display the results for 0.1 < \y\ < 0.7. Non-perturbative 
corrections are determined using PYTHIA with tune A. 

(c) D0 Collaboration (2008) Ant = 0.7 fb^^ [M] 

Measurement of cPa/dpTdy for 50 < pT < 600 GeV in six rapidity regions \y\ < 
2.4 using the Run II iterative midpoint cone jet algorithm (D0 version) with 
R = 0.7 in the i?-scheme. In our comparisons we display the results for \y\ < 0.4. 
Non-perturbative corrections are determined using pythia with tune QW |25j . 

3. Proton-anti-proton scattering at ^/s = 1.8 TeV at the Fermilab Tevatron collider 

(a) CDF Collaboration (2001) dnt = Wph'^ p6] 

Measurement of (Pa/dETdrj for 40 < Et < 465 GeV in the pseudorapidity range 
0.1 < \r]\ < 0.7 using the JETCLU iterative cone jet algorithm with R = 0.7 in 
a special recombination scheme described in equations (3-9) in [26]. Corrections 
for the effects due to the underlying event have been estimated using ambient 
energy measured in minimum bias events and applied to the data. Hadronization 
corrections have not been published. 

(b) D0 Collaboration (2001) Ant = 95 pb^^ [27J 

Measurement of d'^a/dE'rdr] for 60 < Erp < 500 GeV in five pseudorapidity re- 
gions \rj\ < 3.0 using the D0 Run I iterative cone jet algorithm with R = 0.7 
in a special recombination scheme, but with additional corrections of jet rj to 
account for the differences to the i^^^-scheme [28] . In our comparisons we display 
the results for \r]\ < 0.5. Corrections for the effects due to the underlying event 
have been estimated using low luminosity minimum bias data [29] and applied to 
the data. Hadronization corrections have not been published. 

4. Proton-anti-proton scattering at y/s = 546, 630 GeV at the Fermilab Tevatron collider 

(a) CDF Collaboration (1993) at ^ = 546 GeV, £int = 8.58 nb'^ [30j 
Measurement of d'^a/dExdr] for 27.5 < Et < 73 GeV in the pseudorapidity range 
0.1 < \r]\ < 0.7 using the JETCLU iterative cone jet algorithm with R = 0.7 in 
a special recombination scheme, described in equations (3-9) in [2 6 J . Corrections 
for effects due to the underlying event have been estimated using the observed 
calorimeter transverse energy at 90° to the jet axis in CDF dijet events and 
applied to the data. Hadronization corrections have not been published. 

(b) D0 Collaboration (2001) y/s = 630 GeV, at £int = 538 ub'^ [31] 
Measurement of d^a/dExdr] for 21 < Et < 196 GeV in the pseudorapidity region 
\r]\ < 0.5 using the D0 Run I iterative cone jet algorithm with R = 0.7 in a special 
recombination scheme (see [31J). Corrections for the effects due to the underlying 
event have been estimated using Monte Carlo simulations of pp interactions and 
applied to the data. Hadronization corrections have not been published. 

5. Proton-proton scattering at y/s = 200 GeV at RHIC 

(a) STAR Collaboration (2006) Ant = 0.30 pb"^ [32] 

Measurement of d'^a/dpTdr] for 5 < pT < 49 GeV in the pseudorapidity 
region 0.2 < \r]\ < 0.8 using the Run II iterative midpoint cone jet algorithm 
(implemented following [15]) with R = 0.7 in the ii^-scheme. Non-perturbative 
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corrections have not been published. Those were provided to us by one of the 
authors [33]. 



6. Deeply inelastic scattering at ^/s = 318 GeV at HERA 

(a) HI Collaboration (2007) dnt = 65.4 pb'^ \M\ 

Measurement of d'^a/dErdQ'^ for 7 < Et < 50 GeV and 150 < < 15000 GeV^ 
at an inelasticity 0.2 < y < 0.7 in the pseudorapidity region — 1 < ryiab < 2.5 using 
the inclusive jet algorithm in the Breit frame with R = 1 in the S-p-scheme. 
In these comparisons we are using the data 150 < < 5000 GeV'^ for which 
contributions from Z exchange are negligible. The data have been corrected to 
the particle level and for QED radiative corrections, but not for the running of 
a. Hadronization corrections have been determined as the average value using 
DJANGOH [35] and rapgap [36] . 

(b) ZEUS Collaboration (2007) d^t = 81.7pb-^ [37] 

Measurement of d'^a/dErdQ'^ for 8 < Et < 100 GeV and > 125 GeV^ with 
I COS7/1I < 0.65 in the pseudorapidity region — 2 < r/Breit < 1-5 using the inclusive 
jet algorithm in the Breit frame with i? = 1 in the S-p-scheme. In these 
comparisons we are using the data 125 < < 5000 GeV^ for which contributions 
from Z exchange are negligible. The data have been corrected to particle level 
and for QED radiative corrections, including the running of a. Hadronization 
corrections have been determined as the average value using Ariadne \38\ [39] 
and LEPTO [40]. 

7. Deeply inelastic scattering at ^/s = 300 GeV at HERA 

(a) HI Collaboration (2001) dnt = 33 pb^^ [H] 

Measurement of d'^a/dErdQ'^ for 7 < Et < 50 GeV and 150 < < sooOGeV^ 
at 0.2 < y < 0.6 in the pseudorapidity region — 1 < r/iab < 2.5 using the inclusive 
kT jet algorithm in the Breit frame with = 1 in the E'j'-scheme. The data have 
been corrected to particle level and for QED radiative corrections, but not for 
the running of a. Hadronization corrections have been determined as the average 
value from herwig [42j, lepto, and Ariadne. 

(b) ZEUS Collaboration (2002) dnt = 10 pb"^ p] 

Measurement of d^u/dETdQ'^ for 8 < Et < 100 GeV and > 125 GeV^ with 
—0.7 < COS7/J < 0.5 in the pseudorapidity region —2 < r^Brcit < 1-8 using the 
inclusive fcy jet algorithm in the Breit frame with R = \ using the i^^-scheme. 
In these comparisons we are using the data 125 < < 5000 GeV^ for which 
contributions from Z exchange are negligible. The data have been corrected 
to particle level and for QED radiative corrections, including the running of 
a. Hadronization corrections have been determined as the average value using 

ARIADNE, LEPTO, and HERWIG. 
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4 Results 



The ratios of data and theory for ah results obtained in hadron-hadron cohisions and in DIS 
at large four-momentum transfer squared Q^, are shown in Fig. [T]as a function of jet pt- 
For DIS, the variable pT corresponds to the jet transverse energy defined in the Breit frame 
of reference. The uncertainty bars indicate the total experimental uncertainties, defined as 
the quadratic sum of statistical and systematic uncertainties. For purposes of visibility, the 
ratios for different center-of-mass energies are multiplied by factors as given in the figures, 
and the results are sorted from top to bottom in the order of increasing center-of-mass 
energy. 

Figure [2] shows the data points from hadron-hadron collisions, this time as a function 
of XT = 2pt/\/^. For central jet production in hadron-hadron collisions, the variable xt 
is closely related to the hadron momentum fractions xi and X2, carried by the partons. In 
exclusive dijet production with rapidities 2/1=2/2 = 0, the relation xt = xi = X2 holds. 

The comparison of inclusive jet cross section data at central (pseudo-) rapidities as a 
function of xt allows to infer the sensitivity of PDF determinations to the measured data. 

In general, we observe good agreement between theory and data. The data from some 
measurements (CDF at y/s = 546 GeV [30] and ATLAS and CMS at ^ = 7TeV [SHUT]) 
have a tendency of being systematically below the theory predictions. However, the CDF 
and ATLAS data have relatively large uncertainties, and it was discussed in Refs. [T71 [20] 
that within experimental and theoretical uncertainties, the CMS data are in agreement with 
theory. Differences in the shapes between theory and data as a function of pT for the jet 
data at ^/s = 1.8 and 1.96 TeV can be explained by correlated experimental uncertainties, 
as studied in recent global PDF analyses [31 [U H] . 

The theory-data comparisons which are presented in Figs. [1] and [2] for MSTW2008 PDFs 
are also repeated for the PDF parametrizations CTIO, NNPDF2.1, and HERAPDF1.5. 
These results are displayed in Figs. [3] -[HI In most cases we observe a very similar agreement 
as for MSTW2008 PDFs. 

Prom Fig. [1] it is visible that the new LHC jet data have started to go beyond the pT 
reach of the Fermilab Tevatron experiments CDF and D0. Fig. [2] shows that the new CMS 
measurement provides the first jet results which probe the PDFs down to x = 5 • 10~^. In 
this x-region, the PDFs are already constrained by data from inclusive DIS experiments. 
Therefore the comparison of LHC jet data with theory predictions in this overlap region 
will help to understand uncertainties not related to PDFs. It is also visible from Fig. [2] that 
currently the Tevatron jet data still have the highest reach in x. 
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inclusive jet production ^ j^stNLO 

J I rittpV/projects. nepforge.org/fastnlo 

in hadron-induced processes 

(X 4000) Vs = 200 GeV PP 

• STAR cone 0.2 < |ri| < 0.8 



Vs = 300 GeV 
Vs = 318GeV 



DIS 



+ oDAo* ZEUS k-r Q^: from 125 to 5000 GeV^ 

-4^(x200) 

\ ..ATo H1 Q^: from 150 to 5000 GeV^ 

T oDAo* ZEUS k-r Q^: from 125 to 5000 GeV^ 
(x80) 



Vs = 546GeV , , t 

o CDF cone 0.1 <|iil< 0.7 ^ } t DD-bSr 

^^ = 630 GeV t t f M l tiHtf^M } (X25) 

• D0 cone 0.0 <|Ti|< 0.5 T t t .T|TTTTtt tt t + v ; 



Vs = 1 .8 TeV 

o CDF cone 0.1 < |r|| < 0.7 

• D0 cone 0.0 < |ti| < 0.5 

Vs= 1.96 TeV 
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• CMS anti-k^- R=0.5 |y| < 0.5 

all pQCD calculations using NLOJET++ with fastNLO: 

a^(M2)=0.118 I MSTW2008 PDFS | n^ = n, = p.^j^j 

NLO plus non-perturbatlve corrections | pp, pp: Incl. threshold corredtlolns {2-loop) 
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Figure 1: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron collisions and in deeply inelastic scattering at different center-of-mass energies. The 
ratios are shown as a function of jet transverse momentum px- The theory results are 
computed for MSTW2008 PDFs. 
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Figure 2: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron colhsions at different center-of-mass energies. The ratios are shown as a function of 
the scahng variable xt = 2pt/\/s. The theory results are computed for MSTW2008 PDFs. 
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all pQCD calculations using NLOJET++ with fastNLO: 
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Figure 3: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron collisions and in deeply inelastic scattering at different center-of-mass energies. The 
ratios are shown as a function of jet transverse momentum px- The theory results are 
computed for CTIO PDFs. 
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Figure 4: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron colhsions at different center-of-mass energies. The ratios are shown as a function of 
the scahng variable xt = 2pr/y^. The theory results are computed for CTIO PDFs. 
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Figure 5: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron collisions and in deeply inelastic scattering at different center-of-mass energies. The 
ratios are shown as a function of jet transverse momentum px- The theory results are 
computed for NNPDF2.1 PDFs. 
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Figure 6: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron colhsions at different center-of-mass energies. The ratios are shown as a function of 
the scahng variable xt = 2pr/y^. The theory results are computed for NNPDF2.1 PDFs. 
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Figure 7: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron collisions and in deeply inelastic scattering at different center-of-mass energies. The 
ratios are shown as a function of jet transverse momentum px- The theory results are 
computed for HERAPDF1.5 PDFs. 
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Figure 8: Ratios of data and theory for inclusive jet cross sections measured in hadron- 
hadron collisions at different center-of-mass energies. The ratios are shown as a function 
of the scaling variable xt = Ipt/^/s. The theory results are computed for HERAPDF1.5 
PDFs. 



17 



